(N 

o 

(N 



O 






•4— > 



T3 
C 
O 
o 



> 

o 

oo 
\o 

en 

o 



X 



Minimal coupling in oscillator models of quantum dissipation 
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The dissipative harmonic oscillator has two representations. In the first representation the central 
oscillator couples with its position to an oscillator bath. In the second one it couples with its 
momentum to the bath. Both representations are related by a unitary transformation. If the 
central oscillator couples with its position and momentum to two independent baths, no such unitary 
transformation exists. We discuss two possible models of this type and their physical relevance. 
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The Hamiltonian of the harmonic oscillator with mass 
m and frequency uq is characteristically symmetric in po- 
sition and momentum. On the other hand, the dissipa- 
tive harmonic oscillator serves as an example of quantum 
dissipation induced by a gauge field. Let us assume local 
U(l) gauge invariance, requiring expectation values to 
be invariant under the gauge transformation of the wave 
function tp(q) — > exp[iA(q)]tp(q) . This is achieved by in- 
troducing a gauge field A(q) which transforms as A(q) — > 
A(q) — d q A(q). The kinematic momentum 



p^p + A(q) 



(1) 



is gauge invariant. The gauge field A might be decom- 
posed into normal modes and be quantized (hereafter 
fi = 1) 



.4 



1 L^ 



A k 



a k e 



ikq 



a[e- lkq 



(2) 



where a k and a k are bosonic creation and annihilation 
operators of modes with frequency Uk- The parameters 
Afe characterize the strength of the coupling of the gauge 
field to matter. The resulting Hamiltonian reads 



H = 



1 
2m 



P + A(q) 



—q + ^u k a ] k a k , (3) 

k 



with the gauge field minimally coupled to the canonical 
momentum. In the form ([3]) the canonical momentum p 
and the gauge potential A appear as dynamical variables 
and both quantities are gauge dependent. 

We consider in the following the long wavelength ap- 
proximation k — > 0. The space dependence disappears 

e lkq ^ 2 anc l 



A = 



X)*7r( a *- a I)- 

, <^k 
k 



(4) 



A unitary transformation U = exp ( i^4 g J can be applied 
on H . This is the so-called polaron or Goppert-Mayer 



transformation [1]. It brings the Hamiltonian into the 
standard form of the dissipative harmonic oscillator 



H' 



P 
2m 



J n ~2 



E 



Uk 



a-k 



X,, 



Uk 



(5) 



This Hamiltonian was studied first in [2| and later by 
many others, see [3j, |4| and references therein [Jala]- For 
an arbitrary operator O we use the short-hand notation 
\0\ 2 — O^O. Now p is the kinematic (gauge invariant) 
momentum, which in this particular gauge coincides with 
the (generally gauge-dependent) canonical momentum, 
and the particle couples via its position to the modes of 
the gauge field. Thus a specific implementation of quan- 
tum dissipation can be derived from a U(l) gauge invari- 
ance principle. However, in the presence of the gauge 
field, the q-p symmetry under exchange of position and 
momentum is broken, the reason being that we assumed 
gauge invariance which is local in position but not in 
momentum. 

Some authors [7], |8[ considered as well a Hamiltonian 
similar to ([5]), where the oscillator bath couples linearly 
to the momentum 



H" = 



A 2 
P 

2m 



u ~2 



k 



u k 



a-k H P 

Uk 



(G) 



Albeit their formal similarity the physics described by 
H" and H' is quite different. We emphasize that H" 
is inequivalent to H (in the sense of not being related 
by a unitary transformation), even with A as in Eq. (HJ) 
although in both Hamiltonians the oscillator bath couples 
to the momentum. For the coupling to an oscillator bath 
as implemented in H" the term "anomalous dissipation" 
was coined Jj. 

Recently a certain class of dissipative quantum systems 
has aroused interest, where the system couples to two 
heat baths through non-commuting variables [9l4l5|. In 
Refs. |3, llfl, LL5| a spin was considered which couples to 
two independent oscillator baths. In Refs. Ill LlJa llJ] the 
dissipative harmonic oscillator with two competing heat 



baths was investigated. The model was described by the 
Hamiltonian 



Hi 



J q -2 



E 



Uk 



a k -I — q 



^ 2 +E 



f; 



*4 



(7) 



which combines both the conventional coupling of fr and 
the anomalous coupling of H". Such a Hamiltonian can 
be found in real physical situations 1J, |l6[ . Here a k and 



b k are bosonic creation operators of two independent heat 
baths (bath A and bath B) . The frequency of the central 
oscillator is given by loq — ^cu p u! q . In this model the q-p 
symmetry under exchange of position and momentum is 
reestablished, when uj p = cu q and the coupling coefficients 
are identical. 

It is tempting to motivate the model J7) by a gauge 
argument as outlined above for the standard dissipative 
harmonic oscillator and leading to the equivalence of ([3]) 
and (JS|). A second "gauge field" B is introduced which 
couples to the position. The most symmetric and most 
canonical form to introduce this field is via the Hamilto- 



H u = 



^(P + A)\^( i + B 



k 



E^ft 



(8) 



where in the long wavelength approximation B is given 
by B = i^2i —(t>i — bj). However, the Hamiltonians Hi 

and Hu are not equivalent in the sense that they can- 
not be transformed into each other through a canonical 
transformation. As a consequence, the claim made in 
Refs. ll|, [l2|, ll6j that Hi and Hu are unitarily equiva- 
lent is inaccurate. The analysis there performed corre- 
sponds to the Hi model, for which a physical realization 
was found [l6(. The subsequent analysis of the Hamil- 
tonian ([7]) is correct and the results in Refs. [ll|, [l2|, |_16| 
remain unaffected. 

One might naively argue that the Hamiltonian Hi can 
be transformed into Hu by the coordinate transformation 
<2n = q\ ~ Bi, pu = pi- Ai, dfcji = a k j + X k qi/uj k , 
and bui = 6;j + HiQi/vi, where the lower indices I and 
II refer to the corresponding Hamiltonians. But, since 
bi t u and a k u, as obtained by this transformation, do not 
commute, this transformation is not canonical. 

To see this even more clearly one can diagonalize both 
Hamiltonians and verify that their spectra are differ- 
ent. We define the characteristic polynomial Xiui^) = 
det (lj — T^ 1 ' 11 )). Here "H^ 11 ) is the matrix which gov- 
erns the time evolution of the vector whose components 
are the creation and annihilation operators appearing in 
the Hamiltonian Hi (defined as model I in the following) 



respectively in Hu (defined as model II in the following). 
The frequencies of the normal modes are given by the 
positive roots of the characteristic equation Xiui 1 ^) = 0- 
Diagonalising "H^ 11 ) yields 



X5 1 (w)=Wo~ w+-*J p (u/) oj + ^J q (uj) , (9) 

L UJ ALU) 

where we have defined the spectral functions 

J,(w) ee 2^|A fc | 2 <5( W -w fc ) 
fc 

J p (w) ee 2^|/ii| 2 ,5( W -^) (10) 

i 

and, as in [12j, the Riemann-Stieltjes integral transform 

/H =^vf°° y ^-zsgnH^/(M). 

Jo w' \u' - uj z ) I 

(11) 
has been introduced. On the other hand, for model I it 

was found (see [12[) 



Xi X M = 0J q -Jq(u>) U p -J p (ui) 



(12) 



Comparison of Eqs. (JSJ and (TT2l shows that the char- 
acteristic polynomials and therefore the spectra of both 
models are in general different. Importantly, they be- 
come identical if J q (co) or J p (ui) vanishes. 

The question of which of the Hamiltonians Hi or Hu 
is physically relevant has to be answered in favor of the 
Hamiltonian Hi as given in Eq. 0. As already men- 
tioned above, one example for a physical realisation of 
Hi was found in [16j : In a Josephson junction the sys- 
tem which is represented by the particle number differ- 
ence and the phase difference across the junction couples 
to two independent environments given by the bosonic 
exitations of quasiparticle tunnelling across the junction 
and by the electromagnetic vacuum modes. It was shown 
that the corresponding Hamiltonian has the form of H\. 

Another example was given recently by Cuccoli and 
coworkers for spin Hamiltonians [14l . Il7 |. Spin opera- 
tors can often be treated as ordinary canonical variables. 
Then an environmental coupling typically involves both 
the coordinates and the momenta in a symmetric fashion, 



17] 



with no privileged role. This has been shown in Ref. 
for the case of the easy-axis XXZ magnet. 

Although so far it seems to have no direct physical 
application, it is still interesting to analyze model II as 
given in Eq. ([5} • In this brief report we therefore calcu- 
late the equilibrium correlation function of Q = q + B in 
model II and compare it with the equilibrium correlation 
function of q in model I. Our goal is to see if and to what 
extent the physics of both models is different. 

The Heisenberg equations for the bath modes can be 
solved. Plugging the results into the Heisenberg equa- 



tions for Q and P = p + A yields 



which could have been expected from the scaling property 



1 



Q(t) = u p P(t)-uj q / ds K q (t - s)Q(s) + F q (t) 



P(t) = -u q Q{t)-UJ p I dsK p (t-s)P(s)+F p (t). 

— OO 

(13) 
The response kernel is defined as 

poo J f \ 

K n (t) = / -?±-L cos(ujt)duj , n = q,p . (14) 



We defined the force operator 

Fq(t) = *^2 Xkdk exp(-iui k t) + H.c , (15) 

with F p (t) defined accordingly. In Fourier space the 
Hcisenberg equations (fi"3|) read 



(16) 



uj„Q + iuj — —JJuj) P = F q (uj) 
j—-^-J p (u) Q ~ uj p P = F p (u) 



--(ii). 



The equilibrium correlation function C qq (t, (5) = 
^({Q(t), Q(0)})/3 is obtained from equation (fl"6"|) after 
Fourier transformation, namely, 

1 r°° 

C qq l) (t,(3) = - / cMxiiH| 2 cos(u;i)coth(/V2) 

(17) 



[7 p Wo(1 + 7?7p)w + 7 9 w 3 



In the following we set io p 



luq. Moreover we 



focus on the case of Ohmic damping and set J q (Lo) = 
ij q uj and J p (uj) = ij p ui, where 7 g > and j p > 0. The 
poles of the suceptibility lie at 



wo 2 V 



lq-lp 



(18) 



It is seen that at the symmetric point j q = 7 P = 7 the 
real part of the poles is just the oscillator frequency u)q. 
In other words, the frequency remains unchanged by the 
coupling to the baths. No transition to overdamped os- 
cillations can occur. We can compare this with the poles 
u4 I} for model I (see Eq. (11) of Ref. [H). We find the 

simple relation lu± = (l+7 g 7 p )uy- . This indicates that 
the main difference between both models is a dilatation 
of the time scale by a factor s = (1 + 7 g 7 P ) in model I as 
compared to model II. It is seen that when either j q = 
or 7 P = the poles of both models are the same. 

Expression (|17p can be compared with the correspond- 
ing one for model I : C ( qq \t,l3) = \{{q(t),q(O)}) . We 
find the remarkable relation 



C^{t/a,P/s) = sWM), 



<IQ 



(19) 



Xii(suj) = -xi(w) 



(20) 



This also means that the equilibrium mean square ( Q 



in model II is enhanced by a factor s 3 compared with 
(q 2 ) g in model I. In the high temperature limit C qq (t, j3) 
becomes simple (see Eq.(22) of Ref. [12]). At the symmet- 
ric point 7 g = 7 P = 7 and at temperature zero (/3 — > 00) 
it is given by 

C qq \t) ex ^cos(w i/s)exp(-w 7|t|/s) , (21) 



and therefore 



C qq I} (t) ex cos(w t)exp(-w 7|£|) 



(22) 



It seems that although the Hamiltonians Hi and Hu are 
not unitarily equivalent they describe up to a time scale 
factor similar physics. 

In summary we compared two models for a dissipative 
quantum system where a central oscillator couples with 
position and momentum to two independent heat baths. 
One model (model I) is physically well motivated and 
has found by now applications in e.g. Josephson junc- 
tions and magnetic systems. It has been studied before 
in detail 111, |l2j, LLJ] . The other model (model II) has by 



now no obvious physical application and its study seems 
to be a merely academic exercise. A more detailed study 
of this model will be in order if a realistic application for 
it can be found. 

We thank F. Wegner and F. Guinea for useful discus- 
sions. This work has been supported by MICINN (Spain) 
through grants FIS2007-65723 and FIS2010-21372. HK 
acknowledges support from Deutsche Forschungsgemein- 
schaft by the grants K03538/1-2 and from CSIC (Spain) 
through Program JAE. 



* Electronic address: hkohler@icmm.csic.es ; Electronic ad- 
dress: f.sols@fis.ucm.es 

[1] C. Cohen- Tannoudji, J. Dupont-Roc, and G. Grynberg, 
Introduction to Quantum Electrodynamics (John Wiley 
and Sons, New York, 1989), 1st ed. 

[2] P. Ullersma, Physica (Utrecht) 32, 27 (1966); ib. 32, 56 
(1966); 32, 74 (1966); 32, 90 (1966). 

[3] H. Grabert, P. Schramm, and G. L. Ingold, Phys. Rep. 
168, 115 (1988). 

[4] U. Weiss, Quantum Dissipative Systems (World Scien- 
tific, Singapore, 1999), 2nd ed. 

[5] It has been noted in Ref. [a that the Hamiltonian © 
can be slightly modified to become fully invariant under 
translations. 

[6] J. Sanchez-Canizares and F. Sols, Physica A 212, 181 
(1994). 



[7] A. J. Leggett, Phys. Rev. B 30, 1208 (1984). 

[8] A. Cuccoli, A. Fubini, V. Tognetti, and R. Vaia, Phys. 

Rev. E 64 (2001) 066124. 
[9] A. H. Castro Neto, E. Novais, L. Borda, G. Zarand, and 
I. Affleck, Phys. Rev. Lett. 91, 096401 (2003). 
[10] E. Novais, A. H. Castro Neto, L. Borda, I. Affleck, and 

G. Zarand, Phys. Rev. B 72, 014417 (2005). 
[11] H. Kohler and F. Sols, Phys. Rev. B 72, 180404 (2005). 
[12] H. Kohler and F. Sols, New J. Phys. 8, 149 (2006). 
[13] D. D. Bhaktavatsala Rao, Phys. Rev. A 76, 042312 



(2007). 
[14] A. Cuccoli, N. Del Sette, and R. Vaia, Phys. Rev. E 81, 

041110 (2010). 
[15] D. D. Bhaktavatsala Rao, H. Kohler, and F. Sols, New 

J. Phys. 10, 115017 (2008). 
[16] H. Kohler, F. Guinea and F. Sols, Ann. Phys. 31, 127 

(2004). 
[17] A. Cuccoli, A. Fubini, V. Tognetti, and R. Vaia (World 

Scientific, Singapore, 2008), p. 500. 



